Effects of the Nearest-Neighbour Coulomb Interactions on the Ground State of the 

Periodic Anderson Model 
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The magnetic and non-magnetic ground states of the pe- 
riodic Anderson model with Coulomb interaction between f- 
electrons on the nearest-neighbour(NN) sites are investigated 
using a variational method, which gives exact calculation of 
the expectation values in the limit of infinite dimensions. It 
is shown that for a critical value of NN Coulomb interactions 
the magnetic ground state of the periodic Anderson model 
in the Kondo regime is unstable. Factors in terms of the 
physical processes responsible for instability of the magnetic 
ground state are also discussed. Our study indicates the im- 
portance of the NN Coulomb interactions for correlated two 
band models. 

PACS numbers: 71.27.+a, 75.20.Hr, 75.30. Mb 



Over the past decade lot of effort has been devoted to 
the theoretical understanding of the ground state prop- 
erties of the heavy-fermion systems. One of the intrigu- 
ing experimentally observed phenomena in the heavy- 
fermion materials is the variety of magnetic and noft- 
magnetic ground states observed in these materialsaH. 
Most of the theoretical investigations of the magnetic 
properties are done on the basis of the periodic An- 
derson model(PAM) assuming that this model contains 
the essential physics of these materials. Theoretical ap- 
proaches based on the slave-boson techniqueaJH are bi- 
ased towards a paramagnetic ground state while varia- 
tional approaches based on the Gutzwiller method are 
biased towards a magnetic state. These two approaches 
are equivalent in the limit, of large orbital degeneracy. 
Recently, Reynolds et. al£3, studied the magnetic prop- 
erties of the orbitally non-degenerate periodic Anderson 
model using Kotliar and Ruckenstein slave-boson(KRSB) 
formulation of the Gutzwiller method. In this approach 
the Gutzwiller approximation is reproduced at the saddle 
point for T = 0. They found that a magnetic instabil- 
ity exists in the entire Kondo regime and therefore, the 
Gutzwiller approximation is too biased towards the mag- 
netic ground state. The experimental evidence points to 
the gross inadequacy of the existing approaches to de- 
scribe the magnetic behaviour of heavy fermions. 

In addition to the on-site Coulomb interaction in the /- 
band, the other most important interactions which may 
affect the stability of the magnetic ground state of the 



PAM are the on-site Coulomb interaction in conduc- 
tion band and the NN Coulomb interaction in the /- 
band. The influence of the on-site Coulomb interaction 
in the conduction band was recently considered by Itai 
and Fazekasty using the Gutzwiller variational method. 
They found that this interaction reduces the Kondo scale. 
The reduced Kondo scale implies that the transitions 
of electrons from the /-band to the conduction band 
and vice versa, are further restricted by the presence of 
the Coulomb interactions in the conduction band. This 
would lead to further enhancement of the magnetic order- 
ing of the ground state of the periodic Anderson model. 
Consequently, the ground state of the periodic Ander- 
son model with the on-site Coulomb interaction in the 
conduction band would be magnetic in the entire Kondo 
regime. In the presence of the NN Coulomb interaction in 
the /-band all the configurations, having electrons on the 
NN sites are energetically unfavourable and the follow- 
ing physical processes would be operatin (i) /-electrons 
can avoid NN Coulomb repulsion by occupying next to 
nearest neighbour sites. This process is expected to be 
important only when sufficient number of vacant sites 
are available, (ii) Electrons from the /-band may go to 
the Fermi level whereby take advantage of the hybridiza- 
tion interaction to delocalize. (iii) Spin-flip process in 
the /-band through hybridization interaction would also 
lead to energy gain. All these processes would affect the 
magnetic ordering of the ground state of the PAM. The 
purpose of this paper is to investigate the influence of the 
NN Coulomb interaction in the /-band, on the magnetic 
instability of the ground state of the PAM in the Kondo 
regime. To the best of our knowledge, this is the first 
study of the influence of the NN Coulomb interaction on 
the ground state properties of the PAM. 

We consider the extended periodic Anderson model 
given by 

H = J2 e k d L d k« + Y, E f fi f^ + v ^yLf*» + h -°) + 
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where n/j CT = fi a fia, i and j are site indices and k are 
the wave vectors. The first four terms constitute the 



standard PAM and the last term in the Hamiltonian cor- 
responds to the Coulomb interaction between /-electrons 
on the NN sites. Yluj) m ^ e l as * term denotes that the 
sum is taken over NN sites only. The total density of 
electrons n — (J2i<r n / icr + n di<y) /N, where N is the to- 
tal number of lattice sites, is taken to be 1 < n < 2, so 
that there are enough electrons to fill atleast the /-levels, 
and the d-band filling is variable up to half filling. 

To study the magnetic ground state we generalize the 
variational method previously used, £a investigate the 
paramagnetic regime of the PAMifiiil. The general- 
izations are carried out by distinguishing the up and 
down spin electrons in the variational wave function. In 
the previous treatment for the paramagnetic regime of 
the PAM only the lower two spin-degenerate hybridized 
quasiparticle bands were considered, however, for the 
more general case of magnetism, it is required to take 
into consideration all four hybridized quasiparticle bands 
in the variational wave function. 

To investigate magnetic ground state, we choose the 
variational wave function as: 



i^c)=n p *^« 



(2) 



Wk,k>,a,v> l L u la>\°) is the uncorrected 

f[ denotes the product over all occu- 

create quasiparticles in the 



Where \ip uo ) = 

wave function. 

pied states. u ka and l ka 

upper and lower hybridized bands respectively. ll 
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a ^ a 'ka ~ Pkaftcr and u\ a = a ka fl a + Pkad\ a . a ka and 
[3ka are variational functions, which denote the proba- 
bility amplitude for conduction(d) and /-electrons in the 
various quasiparticle bands. The quasiparticle creation 
operators l ka and uj^obey fermion commutation rule if 
a k(T + (3 ka = 1 . The variational functions a ka and (3 klJ 
differ from the choice which diagonalizes the Hamilto- 
nian(Eq. 1) in the absence of Coulomb interactions. Be- 
cause the Coulomb interactions between the / -electrons 
can renormalize the hybridization interaction between d 
and /-electron and thereby can also change the proba- 
bility amplitudes. The correlation operator, Pj, is intro- 
duced to suppress those configurations in the uncorre- 
cted state which are not energetically favourable in the 
presence,of Coulomb interactions. The correlation oper- 
atortL3~t-3 is given by 

Pi = 1 + ^ SvUfia - (1 - d) + ^ S a 1 UfilUfil. (3) 
a a 

The ground state energy per site of the trial wave func- 
tion(Eq. 2) is given by E g /N = (i/> e \H/N\il> e )/{il> e \il> c ). 
The exact calculation of the ground state energy of the 
correlated wave function is not possible since the expecta- 
tion values involve infinite product of operators and one 
needs to adopt some approximate, scheme In this paper 
we use the one-site approximationE3~Ej for calculation of 
various matrix elements appearing in the ground state 



energy per site of \ipc)- The expectation values appear- 
ing in the ground state energy per site of ip c ) , typically 
involve expectation values of the type {..■nfi a nfj a i..) uc . 
In the one-site approximation, such expectation values 
are approximated by 

\...Tlfi ( j'llfj(j>..l uc \...) uc \Tlfio-/uc v^fjcri uc \"')uc \J*J 

where (....) uc — (tJj U c\----\4'uc)- Such an approximation 
implies the collapse of all intersite diagrams in the posi- 
tion space [see Figure |l| . 
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FIG. 1. The collapse of intersite diagrams in the one-site 
approximation. 

The one-site approximation is expected to give exact cal- 
culation of the expectation values in the limit of infinite 
dimensions, since as dimension increases the contribu- 
tion of the intersite diagrams decreases and, vanishes al- 
together in the limit of infinite dimensionsEj. 

Using the one-site approximation described above to 
calculate the expectation values appearing in the ground 
state energy of \ip c ) and minimizing the energy functional 
with respect to the variational functions afc-f, Qfcj_, [3k\, 
Pki by imposing the constraint a ka + (3 ka — 1, the mini- 
mum of the ground-state energy per site is given by 
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where £ fe(T describes four hybridized quasiparticle bands 






(e k +E fa )± (e k ~E frT ) 2 +4V^ 



and {u\ a u k<J ) uc and (lljka)uc correspond to the average 
occupation of the upper {£, ka ) and lower (£, kcr ) quasipar- 
ticle bands. Ef a = Ef — \i a is the renormalized /-level 
energy with the /-electron self energy \x a given by 
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Va = VR a is the renormalized hybridization interaction 
and R a is the renormalization factor. R a and the average 



double occupancy, D of the ground state \ip c ) are given 
by 



R a 



(1-n/) 



.4 



(1 - n/_ (T )(l + s CT ) + dn f - a (l + s_ CT ) 



D = d 2 (l-ri / )n /t n /1 /A 

7 CT = n fa [(l + S l + n f _ a (d 2 -(l + S<J ) 2 )](l - n f )/A 
with A = (1 — n/) + (1 — d 2 )nf^nfi- The density of /- 
electrons ny CT is given by 
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The weight factors $\ and a| ff for /-electrons with spin 
a in the lower and upper quasiparticle bands respectively 
are given by 



aka = 



{e k +E f(T ) + (e k -E frT ) 2 +4V< 



V2 



{e k -E fa ) 2 +W 2 



Pka = 



(e k +E fa ) + (e k -E f 
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The minimization of the ground state energy with respect 
to d yield the following implicit equation for d. 



dd 



n2^ v ° dd 



k a 1 



\ U ka' u ka')uc \^ka' ^ktr'/u 



(e k -E fa ,) 2 +W 2 



-G 



dd 



(9) 



At zero temperatures, we can replace the distribution 
function for the lower and upper quasiparticle bands 
by unit step functions; (l ka l ka )uc = 0( _ £fc CT + v ) an( ^ 
{u\ a u ku ) uc = ®(-£ kcr + v). Here is the unit step 
function, and v is the Fermi level, v is determined by 
fixing the density of total number of electrons per site, 
ii = ^2 a n a . At zero temperatures n a is given by the 
following expression 



= ^E[ e (-^ + -) + e(-^ + -)]. 



(10) 



Before embarking on the numerical calculations it 
would be instructive to compare our approach for the pe- 
riodic Anderson mode]_svith the KRSB reformulation of 
the Gutzwiller methocO. We note that the ground-state 
energy of our variational wave function in the one-site 
approximation and the ground-state energy derived from 
the KRSB approach have different expressions for the ef- 
fective hybridization interaction and the average double 



occupancy of the ground state. The average double occu- 
pancy in the one site approximation and the Gutzwiller 
approximation are given by D and d g (say) respectively. 
If we scale D — > d g in the expression for the effective 
hybridization(V r CT ) in our approach, we find that it re- 
duces to the corresponding expression for the effective 
hybridization in the KRSB method. This further implies 
that the /-electron self-energy (/v), and the average oc- 
cupation of the /-orbitals in both the approaches also 
become the same; thereby the KRSB ground-state en- 
ergy functional and the one-site ground state energy func- 
tional are the same under the scaling of average double 
occupancy of the ground state. Furthermore, since both 
the approaches search for the minimum of the ground- 
state energy in the same physical parameter space, they 
must give the same results at the point of minimum. 
The equivalence of the two seemingly different variational 
methods is surprising. To understand this equivalence, 
we reanalyze the Gutzwiller variational wave function. 
The Gutzwiller wave function has a long history, dating 
back to the wack-of Gutzwiller in 1960's_. The Gutzwiller 
wave fmytionEjE-J is given by \ip g ) = g D \ipo) ■ Recently, 
GebhardllSlij showed that it is more convenient to work 
with the following form for the Gutzwiller wave function. 



IVW = s^o> 



(ii) 



where \<f> ) is an arbitrary normalized one-particle prod- 
uct wave function and K = D — ^\ fj-i^nfia-, where 
/ii CT are the explicit functions of g and the local occu- 
pation of /-orbitals, n/ i(T = (0o|n/io-|0 o ). \ip a ) and \<f> ) 

are connected by \ipo) = ff^ i "' ll<T ™ / "'|^ l o) For the mag- 
netic case the correlator g K can be written as FJ i Qi with 

Qi = 1 + XUfi^Ufii - J2ay<7 n ficn x an d V<? are vari- 
ational parameters which depend on the average occu- 
pation of the /-orbitals. With the redefinition of the 
parameters x and y a , the correlation operators Qi and 
Pi (Eq. 3) are the same. Therefore, the Gutzwiller- 
Gebhard correlator g K in the Gutzwiller approximation 
and our correlator FJ^ Pi in one-site approximation, de- 
scribe the same physics. It is interesting to note that 
the Gutzwiller approximation gives the exact calculation 
of the matrix elements in the limit of infinite dimensions 
and gives identical results as obtained by the one-site ap- 
proximation. Obviously, one-site approximation is much 
more physically transparent and operationally simpler 
than the Gutzwiller approximation. 

Although our variational formalism is valid for arbi- 
trary dimension and dispersion of conduction electron 
band, for simplicity we assume a conduction band with 
a constant density of states p(e k ) = 1/2W lying in the 
energy interval — W < e k < W . 2W is the conduction 
electron band width. We have also taken infinite-?/ limit, 
since at U — oo the ground state of the PAM is strongly 
magnetic with maximum value of total magnetization^]. 
This is an ideal limit to investigate the instability of 
the magnetic ground state in the presence of the near- 



est neighbour Coulomb interactions. In our formalism 
this limit is affected by putting d = throughout i.e., 
by projecting out all the doubly occupied sites. The to- 
tal magnetization, in = ^ ct an a for different values of 
nearest neighbour interaction G, the bare hybridization, 
V and the total electron density, n, is calculated numer- 
ically by solving Eq. (|), Eq. @), and Eq. @ self- 
consistently for n^ fx^, flf\, fi/j,, and v. The numerical 
solution of the self-consistent equations have more than 
one solutions corresponding to strong magnetism, weak 
magnetism and paramagnetism. The relevant solution is 
one with the lowest ground state energy. In the numerical 
calculations we took the conduction electron bandwidth, 
2W = 20eV and the /-level, E f = -l.beV below the 
middle of the conduction band. We have taken this par- 
ticular choice of parameter values for reasons of compar- 
ison with earlier work of Reynolds et all3 in the absence 
of nearest neighbour Coulomb interaction between the 
/-electrons. 

In Figure 2, we have plotted the total magnetization as 
a function of the NN Coulomb interaction G, for the total 
density of electrons, n=1.95 and n=1.9. For G/\Ef\ = 0, 
the ground state is strongly ferromagnetic with total 
magnetization, m = 0.96 for n = 1.95 and m = .94 
for n = 1.9. With increasing value of G the magneti- 
zation decreases upto to a critical value of G, where we 
see a crossover from strong ferromagnetism to weak fer- 
romagnetism with total magnetization m = 2 — n and 
then from weak ferromagnetism to paramagnetism. Fig- 
ure 3 shows the magnetic phase diagram of the extended 
periodic Anderson model. 

In Figure 4 we have plotted the hybridized quasipar- 
ticle bands ^ for strongly ferromagnetic {G/\Ej\ = 0), 
weakly ferromagnetic (G/\Ef\ = 2) and paramagnetic 
(G/\Ef\ = 5.6) ground states. We find that due to the 
renormalization of the hybridization interaction and the 
/-electron energy there is a redistribution of the den- 
sity of states and to accommodate the redistribution of 
density of states the Fermi level also moves to keep the 
total density of electrons fixed. For the strongly mag- 
netic ground state the Fermi level lies in the lower down 
spin hybridized band (Cfei) an d the upper up-spin hy- 
bridized band(£i.). The lower up-spin hybridized band 
(£JmO is completely full. At a critical value of G/\Ef\ 
all the electrons in the upper hybridized up-spin band 
are transferred to the lower hybridized down-spin band. 
Then we see a crossover from strongly ferromagnetic to 
weakly ferromagnetic ground state. In the weakly fer- 
romagnetic ground state, the lower hybridized up-spin 
band is completely full with total density of up-spin 
electrons, n-\ = 1, therefore the total magnetization, 
m = n^ — (n — n-f) = 2 — n. The total magnetization re- 
mains unchanged in the entire weak ferromagnetic regime 
till the fermi level also lies in the lower hybridized up- 
spin band. To understand the magnitude of jump in the 
magnetization at the point of crossover from strong fer- 
romagnetism to weak ferromagnetism and then another 



crossover from weak ferromagnetism to paramagnetism, 
it would be instructive to calculate the density of states 
p„{uj) of hybridized bands, £ fc(T . It is given by 



p±H = V<5(, 



k 



&) = £*("-**) 



d& 



de h 



(12) 



where Xi's are the roots of uo — £ feo . = 0. At a given en- 
ergy w, the density of states is proportional to the slope 
of the hybridized bands at u). The crossover from strong 
to weak ferromagnetism is due to shift of the Fermi level 
from the upper hybridized up-spin band with larger den- 
sity of states to the lower hybridized down-spin band with 
smaller density of states, resulting in a large decrease in 
the number of up-spin electrons. In the crossover from 
weak ferromagnetism to the paramagnetism, the change 
of the Fermi level in the lower hybridized up-spin band to 
the lower degenerate up and down spin bands acompanies 
a relatively weaker change of the density of states. 

Figures 5 shows how the NN Coulomb interaction 
renormalizes the average occupation of the /-orbitals. 
In the absence of the NN Coulomb interaction the to- 
tal magnetization, m ~ n/f. In the strongly ferromag- 
netic regime (G/ 1 Ef \ < 2.0), for small values of G, up-spin 
electrons from the /-band are transferred to the conduc- 
tion band through the hybridization interaction. With 
increasing values of G, more and more vacant sites are 
available and spin flip processes through hybridization 
interaction becomes energetically favourable. Therefore, 
we see an increase in the number of down-spin electrons 
and a decrease in the total magnetization. The weak fer- 
romagnetic regime is stabilized by energy gain through 
the transfer of electrons from the /-band to the conduc- 
tion band and by occupying next to nearest neighbour 
/-electron sites, since there are sufficiently large number 
of vacant sites available in this regime. 

In this paper we have investigated the magnetic and 
non-magnetic ground states of the extended periodic An- 
derson model, using a variational method based on the 
one site approximation. We have shown through the cal- 
culation of the magnetic phase diagram that for U = oo, 
the non-magnetic ground state is stabilized above a crit- 
ical value of nearest neighbour Coulomb repulsion be- 
tween the /-electrons. The one site approximation used 
in this paper gives exact calculation of the matrix ele- 
ments in the limit of infinite dimensions. It will be very 
interesting to investigate the magnetic properties of the 
periodic Anderson model by including the (dimension) -1 
contributions through two-site approximationEl Cer- 
tainly, it is desirable to extend our calculations to study 
antiferromagnetic ground states also. 
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FIG. 3. Magnetic Phase diagram of extended •periodic An- 
derson model 
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FIG. 2. The total magnetization as a function of the 
Coulomb repulsion between f -electrons on the nearest neigh- 
bour sites. Here we have taken the on-site Coulomb interac- 
tion to be infinitely large and bare hybridization, V = 1. 
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FIG. 5. The average occupation of the f-orbitals n/ CT are 
plotted as a function of G/\Ef\ for n — 1.95. 



FIG. 4. Hybridized quasiparticle bands (^ ) as a function 
of conduction band energy(ek) for the total density of elec- 
trons, n = 1.95 and V — 1.0. The j and J, correspond to 
the hybridized up-spin band and hybridized down-spin band 
respectively. For G/\Ej\ = 5.6, the hybridized band are spin 
degenerate. 



